Supercurrent-induced temperature gradient across a nonequilibrium SNS Josephson 

junction 
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Using tunneling spectroscopy, we have measured the local electron energy distribution function 
in the normal part of a superconductor-normal metal-superconductor (SNS) Josephson junction 
containing an extra lead to a normal reservoir. In the presence of simultaneous supercurrent and 
injected quasiparticle current, the distribution function exhibits a sharp feature at very low energy. 
The feature is odd in energy, and odd under reversal of either the supercurrent or the quasiparticle 
current direction. The feature represents an effective temperature gradient across the SNS Josephson 
junction that is controllable by the supercurrent. 

PACS numbers: 74.50.+r, 73.23.-b, 85. 25. Am, 85.25.Cp 



The study of nonequilibrium phenomena in supercon- 
ductors has a rich history extending over several decades 
0. Advances in microfabrication techniques in the 
1990's have broadened the range of possible experiments. 
In 1999, Baselmans et al. [2( turned a superconducting- 
normal-superconducting Josephson junction (SNS JJ) 
into a "7r-junction" by driving the electron energy dis- 
tribution far from equilibrium with a normal control cur- 
rent applied perpendicular to the supercurrent in a four- 
terminal geometry. Later, Huang et al. 3] produced 
similar behavior in a three-terminal sample with the con- 
trol current injected from a single normal reservoir. The 
physical explanation behind the 7r-junction 4] 5] is the 
same in the two geometries: the electron energy distri- 
bution in thepresence of the control current has a double- 
step shape which selectively populates the Andreev 
bound states that carry current in the opposite direction 
to that in equilibrium. 

While the physics of the three-terminal and four- 
terminal 7r-junctions is similar, the current flow and en- 
ergy distributions are different in the two geometries. In 
the four-terminal experiment, the control current and su- 
percurrent coexist only in the region where the two wires 
cross. In the three-terminal experiment, the wire joining 
the two superconductors carries both supercurrent and 
normal (quasiparticle) current Q. It was predicted in 
Ref. |3 that this coexistence of quasiparticle current and 
supercurrent will lead to a left-right asymmetry in the 
effective temperature, due to mixing of the even and odd 
components of the distribution function, f(E). In this 
process, supercurrent redistributes the Joule heat due 
to the quasiparticle current between the different parts 
of the system. As a result, the double-step distribution 
function is distorted, and in the absence of inelastic scat- 
tering, this distortion is approximately proportional to 
the energy spectrum of the supercurrent-carrying states. 
We present here the first measurement of the distortion 
in f(E), using a local tunnel probe. 

Fig. la shows the three-terminal sample geometry. A 
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FIG. 1: (a) SEM micrograph of a sample and schematic mea- 
surement circuit diagram, (b) Spectral supercurrent js for 
4> = 7r/2 calculated numerically for the sample shown in (a), 
(c) dl/dV of the tunnel probe (P) measured at B=0.4 T, 
where superconductivity is suppressed. The 10% reduction in 
dl/dV at low bias is due to dynamical Coulomb blockade, as 
discussed in the text. 



normal metal wire connects two superconducting elec- 
trodes (Si and S2) and a normal reservoir (N). The su- 
perconducting electrodes are kept at the same chemical 
potential, which we define as the zero of energy. With the 
normal electrode biased at a potential Un, the distribu- 
tion function f(E) can be calculated easily in the absence 
of supercurrent, energy relaxation, and proximity effects. 
In the normal reservoir, f(E) — fFD{E+eUisf) is a Fermi- 
Dirac function displaced in energy by —eUjf. Define the 
even and odd (in energy) components of f(E) as follows: 
ME) = 1 - f(E) - f(-E) and f L = f(—E) - f(E). 
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The boundary conditions at the NS interfaces for ener- 
gies below the superconducting gap, A, are f T = and 
df]j/dx — 0, assuming high-transparency interfaces, no 
charge imbalance in the superconductors, and no heat 
transport into the superconductors |9j. The solution 
for f(E) at the NS interface is the double-step shape: 
f(E) = 0.5 * [f FD (E + eU N ) + f FD (E - eU N )}. 

In the presence of proximity effects and supercurrent, 
the kinetic equations become more complicated [Tof : 



— = 0, jT = D T (x)—+j E f L +T( X )—; (la) 



-=0,^ = ^)^+^-11%; (ib) 

The spectral charge and energy currents, jr(E) and 
jjj(E), are conserved only in the absence of energy re- 
laxation. The energy-dependent coefficients Dt, Dli .7 en 
and T can be calculated from the Usadel equation 
and vary with the superconducting phase difference <j) 
between Si and S 2 - Fig. lb shows the spectral supercur- 
rent density, j E , for <f> = tt/2 obtained from a numerical 
solution of the Usadel equation appropriate for our sam- 
ple. In the presence of the proximity effect, the kinetic 
equations must be solved numerically. But one can un- 
derstand the origin of the effective temperature gradient 
predicted in Ref. 8 by noting that the terms with je mix 
/t and fi, in the kinetic equations. If we neglect the 
T terms and ignore the energy dependence in Dt and 
Dl, then Eqs. (la) and (lb) are easily solved analyti- 
cally. The solution shows that Jl along the horizontal 
wire connecting the two superconductors contains a spa- 
tially antisymmetric contribution proportional to je- 

The experiments described in this paper were per- 
formed on four samples, with similar results obtained on 
all four. Most of the data reported here were taken from 
the sample shown in Fig. la, with Ri = R 2 = 7.0 f2, 
R3 = 16.9 ft and tunnel resistance Rt = 115 kil. The 
samples were fabricated using a PMMA stencil patterned 
by electron-beam lithography and multiple-angle evapo- 
rations [l^ . The three-legged normal Ag wire is 70 nm 
wide and 30 nm thick. The superconducting Al elec- 
trodes are 60 nm thick, separated by a distance of 1.4 /Ltm. 
The normal reservoir is 200 nm thick. The tunnel probe 
(connected to the pad labelled P) is 20 nm-thick Al oxi- 
dized in 80 mbar of Ar(90%)-O 2 (10%) for 6 minutes be- 
fore deposition of the Ag wire. 

We use tunneling spectroscopy to measure the shape of 
the distribution function of the Ag wire at a position close 
to the superconducting reservoir Si, in a manner similar 
to that in Ref. 0. A complication in our experiment is 
that the density of states in the Ag wire is strongly mod- 
ified by proximity to Si and S 2 - The current-voltage 



characteristic of the probe tunnel junction is 

I(V) = - JL y dEn s (E) J dsP(e) (2) 

[fs(E)n N (E -eV-e)(l- f N (E - eV - e)) 
- (1 - fs{E))n N {E -eV + e)f N (E - eV + e)] 

where Rj> is the normal state tunnel resistance, and 
ns are the normalized densities of states, and /n and 
fs are the electron energy distribution functions on the 
Ag and Al sides of the tunnel junction, respectively. Eq. 
(J2J includes the effects of "dynamical Coulomb block- 
ade" ^3 through the quantity P(e), the probability for 
an electron to lose energy e to the resistive environment 
while tunneling through the barrier. 

Our strategy was to perform a series of experiments in 
which all of the quantities in Eq. (J2J) are known except 
for one, and then deconvolve the data to find the missing 
quantity. The first step is to find P{s), which can be de- 
termined from dl / dV (V) data with a large magnetic field 
B to suppress superconductivity. Fig. lc shows dl/dV 
vs. V with B = 0.4 T. In agreement with previous work 
[l4T |. dl/dV depends logarithmically on V over nearly 
three decades of voltage. For a tunnel junction in a resis- 
tive environment R, with R <C Rk = h/e 2 , the function 
P(e) is well approximated by: 



P(e) 



Eq yEoj 





,e < E Q 
,£ > E 



(3) 



where a = 2R/R K and E = e 2 /{iraC). The dl/dV data 
in Fig. lc can be fit to this form for P(e) with a — 0.0176 
and Eq = 2.0 meV, which correspond to the reasonable 
values R = 230 fl and C = 1.4 fF. 

The second step is to determine niq{E) 1 the density 
of states of the Ag wire in close proximity to Si. In 
equilibrium, fs(E) and ffj(E) are Fermi-Dirac functions 
at a temperature close to that of the cryostat (35 mK). 
In most of our samples, the dI/dV(V) data in zero ap- 
plied magnetic field show an anomalous sharp feature 
at eV w A, which disappears when B > 10 mT (lij: 
hence, all the low-field data reported in this paper were 
taken with B = 12.5 mT. We use the standard BCS 
form for ns(E) of the Al probe, with a small depairing 
parameter to account for the applied magnetic field [16j . 
To determine tin{E) without ambiguity, we must know 
the gap A of the Al probe. For that purpose we fabri- 
cated, on the same chip, a reference SIN tunnel junction 
with no proximity effect on the N side. The dl/dV data 
for the SIN junction, shown in the right inset of Fig. |2 
could be fit perfectly with Eq. (0) assuming a BCS form 
for ns(E) with A = 274 fieV, a depairing parameter of 
7 = T/A = 0.0020 and with Dynamical Coulomb Block- 
ade parameters determined from high-field measurements 
on that junction. The rather large value of A is common 
for thin thermally-evaporated Al films 01- dl/dV data 
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for the sample tunnel junction, shown in the left inset 
to Fig. 2, looks noticeably different from the dl/dV data 
on the reference SIN tunnel junction due to proximity 
effect in the Ag wire. With A determined, the dl/dV 
data can be deconvolved to produce un{E) of the Ag 
wire, shown in the main panel of Fig. 2. Equilibrium 
dl/dV measurements were repeated with supercurrent 
passing through the SNS J J, since Un(E) depends on <j>, 
the phase difference between Si and 52- Fig. 2 also shows 
iin(E) obtained with Is = 0.97 c and — 0.9J C , which are 
indistinguishable from each other. 
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FIG. 2: (Color online) The n N {E) of the Ag wire at the 
position of the probe electrode, with (red) and without (black) 
supercurrent passing through the SNS Josephson junction, 
obtained from dl/dV data such as shown in the left inset. 
The solid line is a fit to the solution of the Usadel equation 
discussed in the text. The (blue) dashed line shows un{E) 
needed to fit nonequilibrium data with Un = 53 fiV. Right 
inset: dl/dV data on a reference SIN tunnel junction, used to 
determine the gap A of the superconducting probe. 

To obtain a decent fit of the data in Fig.|3to the den- 
sity of states calculated from the Usadel theory, we need 
to assume a considerably smaller value of A ks 150 /xeV 
in the superconducting reservoirs feeding the supercur- 
rent than that measured with the SIN reference junc- 
tion. The reason for this smaller value is the presence 
of a normal-metal-superconductor bilayer close to these 
reservoirs (see Fig. 0. Once we know this A, we find 
the Thouless energy Et = HD/L S characterizing the 
normal-metal wire of length L$ between the supercon- 
ductors by fitting the temperature dependent supercur- 
rent at Un = analogously to Ref. 3. The quality of the 
fit is fairly insensitive to the precise value of A, and we 
obtain Et = 5.56 fieV. In the fit, we also corrected for 
the finite size of the superconducting terminals and cal- 
culated the position-dependent order parameter A self- 
consistently close to the NS interface. Once we know 
these two energy scales, we can calculate the local den- 
sity of states un(E) at the position of the probe. As 



there is some ambiguity about the position due to the 
finite normal-metal-superconductor overlap, we fit also 
this within the limits set by the sample geometry. This 
fit is shown as a solid line in Fig. [21 

The third step in the experiment is to apply a voltage 
Un to the normal reservoir, and repeat the dl/dV mea- 
surements across the probe junction with various values 
of Un and Is- Fig. [3] shows three sets of dl/dV vs. V 
data with Un — 53 /iV, and with Is = 0, — 0.9J C , and 
+0.9I C . The double peaks in dl/dV near V = ±A arise 
from the double-step shape in /n(E), the distribution 
function in the Ag wire. Since we know P(e), ns(E), 
and un(E), and since fs(E) is a Fermi-Dirac function, 
we planned to determine Jn (E) from a deconvolution of 
these dl/dV data. Instead, we found that we were un- 
able to fit the data shown in Fig. [3] using Eq. , with 
any physically permissible form for /n{E). Only by al- 
lowing hn{E) to vary slightly from the form obtained 
in equilibrium (dashed line in Fig. 2), could wc fit the 
nonequilibrium data. 
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FIG. 3: dl/dV data of the probe tunnel junction with the 
normal reservoir biased at Un — 53 /iV. Inset: Enlarged view 
of the double-peak at positive bias, showing the differences 
between the data sets taken with Is = and ±0.97 c . 

Fortunately, we can extract the supercurrent- 
dependent changes in Jn{E) by considering differences 
between data sets with opposite values of Is- The 
equilibrium measurements shown in Fig. 2 confirm our 
expectation that un(E,Un — 0, Is) = un{E,Un — 
0, —Is)- This symmetry holds approximately also 
for Un i= 0; hence the quantity dI/dV(V,UN,Is) — 
dI/dV(V,UN,—Is) depends only on the difference 
Sf N (E) ee f N (E,U N ,I s ) - f N (E,U N ,-I S ) = Sf L (E), 
the change in the function Jl(E). The left panel of 
Fig. 01 shows 5/n defined in this way, with Un = 22 /xV 
and Is = 0.97 c . There are sharp peaks at E w ±15// V, 
odd in energy. This feature in Jn{E) represents changes 
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in the local effective electron temperature, according to 
Eq. (5) in Ref. To confirm that the peaks are not 
due to a geometrical asymmetry in the sample, such as 
a flaw in the fabrication process, a second determination 
of S/n(E) is shown superposed on the first. This second 
form was obtained by reversing the signs of both Un and 
Is: ME, -U N , -I s ) - ME, ~U N , I s ). The two data 
sets are in near perfect agreement. 
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FIG. 4: (a) Sf{E) = f N (E,U N ,I s ) - f N (E,U N ,-I s ) for 
Un = 22 fiV and 7s = 0.97 c . (b) Same quantity for Un = 
53 /TV and 7s = 0.97 c , with the SNS J J in the 7r state. Notice 
the sign change between the two data sets. In both figures, 
a second data set is shown with the signs of both Un and 7s 
reversed. Solid lines are numerical solutions to Eq. (1). 

If we increase Un, we can convert the SNS J J to a ir- 
junction 0. The SNS J J crosses over from the 0-state 
to the 7r-state at Un — 30 fiV. The right panel of Fig. 0] 
shows S/n(E) data in the 7r-state for Un = 53 fj,V and 
7s = 0.97 c . Compared to Fig. 4a, the sign of the low- 
energy feature in SfN{E) is reversed, demonstrating that 
the phase difference <j), rather than the supercurrent, de- 
termines the sign of the effective temperature gradient 
across the SNS J J. 

Let us compare these results to those obtained from 
Eqs. (JIJ using the parameters Et, A obtained from the 
previous fits, with no additional fitting. The computed 
theory curves are shown as solid lines in Fig. 01 and they 
agree with the experimental data rather well. Note that 
for this calculation, we entirely neglected the effect of 
inelastic scattering. 

As stated above, we have recovered 8/n(E) without 
knowing the exact form of Hn(E). We have verified 
that our results are rather insensitive to small changes 
in un{E) by performing the deconvolution using several 
different forms for un(E), with the value of iin(0) vary- 
ing from 0.1 to 0.5. The amplitude of the oscillations in 
6fN{E) are inversely proportional to tin{E) near E = 0, 
but the overall shape does not change from those shown 
in Fig. 0] To extract the full Jn (E) from our data would 
require a separate measurement of tin(E) (using a nor- 
mal tunnel probe [18J) at each value of Un and Is used 
in our experiment. As Un increases, deviations of Un(E) 
from its equilibrium form may arise from nonequilibrium 
processes in the superconducting reservoirs, where quasi- 



particle current is converted to supercurrent [l9|. 

The feature we have measured has been called a 
"Peltier-like" effect in Ref. 0, because an applied current 
induces a temperature gradient. It should be emphasized, 
however, that our experiment is far from equilibrium. We 
are not measuring a linear response transport coefficient, 
as was done with the thermopower of Andreev interfer- 
ometers [13, El R is unlikely that one could achieve 
real local cooling in our experiment, as one does with the 
conventional Peltier effect. Perhaps the most intriguing 
aspect of our results is that the nonequilibrium feature in 
f(E) shown in Fig. 4 provides an approximate visualiza- 
tion of the spectral supercurrent density ]e - a quantity 
whose energy dependence has previously been inferred 
only indirectly through 7r-junction experiments. 
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